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The many-body correlation effects in the spatially separated electron and hole layers in the cou- 
pled quantum wells (CQW) are investigated. The specific case of the many-component electron-hole 
system is considered, v S> 1 being the number of components. Keeping the main diagrams in the pa- 
1 rameter 1/v allows us to justify the RPA diagram selection. The ground state of the system is found 

to be the electron-hole liquid which possesses the energy smaller than the exciton phase. The con- 
nection is discussed between the results obtained and the experiments in which the inhomogeneous 
q ' state in the CQW is revealed. 
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I. INTRODUCTION 



^ ■ The investigation of spatially separated electrons and holes in the coupled quantum wells (CQW) is motivated by 
expecting that the electron-hole pairs forming a long-living bound state, namely exciton, can experience the Bose- 
Einstein condensation^. The interest in such systems has greatly grown in the recent years due to the increasing 
ability to manufacture the high quality quantum well structures in which electrons and holes are confined in the 
different regions between which the tunneling can be made negligible 2 . For this reason, in the CQW the exciton 
lifetime becomes by several orders of the magnitude longer than the lifetime of excitons in bulk materials, enabling 
^ \ an experimental observation of the exciton Bose-Einstein condensation. 

As early as decade and a half, Lozovik and Berman, using the variational approach, predicted the existence of an 
electron-hole condensate phase in CQW—. The further theoretical investigations of such systems reveal that the phase 
O ■ diagram of the system can be rather ricbA 

The main goal of the paper is to find out how the many-body correlations in the electron-hole system in the CQW 
influence its ground state. A specific case of the many- component spatially separated electron-hole plasma in the 
' CQW is studied in the current paper. It is assumed that v different kinds of electrons arc confined in one layer 
CNj . of the CQW, while v different kinds of holes are confined in the other layer, the parameter being v 3> 1. For the 
C*~) ■ first time, the large v approximation was proposed for the investigation of a 3D electron-hole liquid in many- valley 
\ semiconductors^. Then, this approach got a further development^. 

In the model under consideration, the thickness of the layers is supposed to be so small that the carriers are the 
2D ones. Let as be the effective Bohr radius and I be the inter-layer distance. Below, we mainly consider the case 
I ^> as- Under such condition, the isolated electron- hole pair may form an exciton possessing the in-plane radius 
Rex ~ clb (1/o.b) 3 ^ 4 3> as 1 ■ Let n be the carrier concentration. The concept of excitons has sense until the average 
in-plane inter-exciton distance being of the order of n" 1 / 2 is larger than R ex , i.e. R%.n <C 1. At higher density the 
exciton system brakes down into the electron-hole one. 

It is shown in this paper that the negative contribution of the energy induced by the intra-layer carrier correlation 
strongly reduces the total energy. This effect facilitates the electron-hole liquid formation. For this phase, the 
equilibrium carrier density n eq , which is the same for the both layers, is found to obey the condition R\ x n eq ~ 1. Thus, 
the criterion for an existence of the exciton state is violated. Therewith, the electron-hole liquid is shown to possess 
the energy lower as compared to the exciton state. Note, that the similar conclusion was made for conventional 3D 
semiconductors, as well££. 

The homogeneous electron-hole liquid state predicted in this paper is realized within the assumption that the 
carriers possess the infinite lifetime. If not, the system breaks down into the low-density exciton gas phase and the 
dense electron-hole liquid phase. Between these phases the dynamic equilibrium is established which is supported by 
the external pumping. A similar scenario was proposed by Keldysh for conventional 3D semiconductors^. 

The diagrammatic technique is used to calculate the correlation contribution to the energy of the electron-hole 
state. The approach is based on the 1/v expansion, what allows one to justify the random phase approximation 
(RPA) diagram selection. It is interesting to note that, for I 3> as, the ground state energy as well as the equilibrium 
carrier concentration n eq do not depend on the parameter v. This unexpected feature allows us to expect that the 
results obtained may be, at least qualitatively, extended even to the region v ~ 1. 
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The paper is organized as follows. First, we describe the model of the electron-hole system in the CQW. Then, the 
self-consistent diagrammatic approach is proposed to calculate the electron and hole Matsubara Green functions and 
estimate the chemical potential /j, for the electrons and holes at temperature T <C ef, £f being the Fermi energy for 
both electrons and holes. This allows us to find the equilibrium electron-hole liquid state concentration n eq . Though, 
to obtain the main results, it is assumed that / ^> attire case I < as is also analyzed at the end of the paper. In 
conclusion, we consider a possible relation between the results obtained and the experiments, in which the nonuniform 
state of the system in CQW is revealed^—. 



II. MODEL AND THE CALCULATION OF THE SELF-ENERGY 



To describe the spatially separated electron-hole system in the CQW, it is supposed that the electrons are confined 
within one infinitely narrow 2D layer, while the holes are confined within the other layer. For the sake of simplicity, 
it is supposed that the effective mass of both electrons and holes are the same. Below the system of units is used 
with the effective electron (hole) charge e = 1, effective electron (hole) masses m = 1, and the Planck constant H = 1. 
Then, the effective Bohr radius is a# = 1, and the energy is measured in the Hartree units. 

Let n be concentration of electrons or holes. According to the model under consideration, the Fermi momentum 
Pf = 27T 1 / 2 (n/v) 1 / 2 and the Fermi energy ep = 2ixnjv are the same for both electrons and holes. We confine ourselves 
to the case of strongly degenerated plasma, that is, the temperature T <^ef- 

The Hamiltonian of the system is H = Ho + U, Hq being the kinetic energy and U being the Coulomb interaction. 
In the second quantization one has 

H = Y<Y a « Ak)aaAk);d = W ^ ^(Ik|)a+ CT (k)a+ CT ,(k')a a ^(k'-q)a a(T (k + q). (1) 

acrk aa'ercr'kk'q 

Here a — e for the electrons and a = h for the holes. Then, a^ a (k) and a aiy (k) are the creation and annihilation 
operators, k is the in-plane 2D momentum, a = 1, ...v denotes the kind of the electron or the hole component. The 
Coulomb interaction in the momentum representation reads 

Let G° a a (uj, p) = (iui + yt* — s p ) 1 be the electron (hole) free Matsubara Green function and G a . a (u, p) = 

{iuj + [i — S a (7 (cj, p) — Ep)" 1 be the electron (hole) total Green function. Here u) is the Matsubara frequency, 
e p = p 2 /2, S Q!CT (cli,p) is the mass operator, and \i is the chemical potential. Let us consider the system of self- 
consistent equations whose diagram representation is shown in Fig[Tj 




FIG. 1: The self-consistent diagrammatic equations for the Green function and mass operator. Here E is the mass operator. 
The double solid line denotes the total Green function G. The wavy line represents the interaction (J2]) and the double wavy 
line corresponds to the renormalized interaction. 



The first diagrammatic equation in FigJT] reads 

Ea,a(e,p) = EW (e,p) + E« (e,p). 



E& (e,p) = % E £W(0)GW(k,u;), 



V 



Ei 2 l(£,p) 



T 
V 



^2 Uaa ( k > w ) G <*° (P + k, e + w) . 



(3) 
(4) 

(5) 
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Here T is the temperature, and V is the square of the layers. 

The renormalized interaction U aa (p, uj) in Eq. ([S]) obeys the second diagrammatic equation in FigfT]and reads 

U aa (p, Uj) = U a a (p) + U aa (p) Uaa (P, £) U aa (p, w) ■ (6) 

The polarization operator 

II (p, e) = U a<7 (p, e) = -^] Gaa (p + k,e + uj) G aa (k, w) (7) 

does not depend on the kind of the particles. 

Let us explain the rule of diagram selection. Each fermion loop contributes a factor v to the diagrams. Among all 
the diagrams of the given order, only the diagrams are retained which contain the maximal number of the fermion 
loops. Thus, our approach is a 1/v expansion which results in keeping only the diagram of the RPA kind. The 
diagrams for which the bubbles are linked by more than one interaction line are small as compared with the RPA 
ones in the parameter 1/V.For the same reason, the diagrams are small, which contain any vertex corrections. 

The expression for the renormalized interaction follows from Eq. ([5]) 

U aa (p,uj) = U(p,uj) = - !! aa ( iP L ( r- (8) 

\-V- Uaa (p)n(p,w) 

Let us note that, like the polarization operator H a a (p,w), both the mass operator T, a ,a (£,p) and the renormalized 
interaction U aa (p, w) do not depend on the kind of the particle. 

Our first goal is to calculate the chemical potential ^ as a function of the density n using the self-consistent system 
of Eqs. (JSI-jH]). Let us remind the exact relation 

fi-S(0,p F )=p F /2. (9) 

Then, the convergence method is used. If the interaction is neglected, the mass operator E (0,pf) ~ and, therefore, 

fj, m fj, = pp/2 = 2im/ v. (10) 

To find the next iteration for the chemical potential, let us calculate E(0,p^) taking /j, — /io- I n this case, one 
should replace G by Go in the above system of equations. First, taking into account Eq. we obtain 

E&(e,p) = ! E Uaa-{Q)G a , a ,{^u) = 2TTln. (11) 

Then, to find the next contribution to the mass operator E^ CT (uj, p) , one should first estimate the renormalized 
interaction U aa (k, uj). For this purpose, let us consider the polarization operator in zero order in the interaction: 

jj, c , 1 f ,2 h n ( k ) ~ n ( k + P) / 10A 

n(p,e) w 7—2 / a k— — — — -2— , (12) 



(2tt) z J ihe + k 2 /2 - (k + p)^ /2 



where no (k) 



exp 



(fc 2 /2- M0 ) + 1 



-1 



is the Fermi distribution function. 



For small momenta and frequencies p < p F , e ^ £f, one finds IIo (p,e) — — 1/2tt. Substituting this value into Eq. 
((HI) gives the estimate 

U (P, s) = — i c 1 « 1. (13) 

1 + Uaa (P) f 5¥ V 

This estimation is valid if 



what is supposed below. 



n x l z lv « 1, 



(14) 
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Let us turn to the behavior of the renormalized interaction U(p,e) for large momenta |p| 3> Pf- First, let us 
estimate the polarization operator (fT2"j) for such momenta. Consider, for example, the contribution II' (p, e) associated 
with the factor n aa (p + q) . It is evident that momenta k alone contribute to the polarization operator for which 
|p + k| < pp. For this reason, |k| 3> pp. Therefore, (p + k) <C k 2 ss p 2 and 

• , . 1 nlv 

IT (p,e) - 



(2tt) 2 ife + p 2 /2' 

Similarly, one can estimate the contribution IIq (p, e) connected with the term n p in Eq. (|12|) . As a result, we obtain 

1 np 2 

(27r) 2 ^e 2 + (p 2 /2) 2 

Taking into account (|15p, one can estimate the renormalized interaction 



n(p,e)w— — 2 — — 2, P>PF- (15) 
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Substituting (p~3|) into ([5]) gives the estimate for the contribution to the self-energy Hal (w, p) coming from |k| < pp. 

^ 2 )( £ ,p)«i4. (17) 
Let us turn to the contribution to S^j (w,p) gained from |k| 3> pf- It is convenient to represent it in the form 

K$(e,P) = ~ E (f/(k)GL(p + k, £ + a,)- 

~A^(k)G^ <r (p + k,e+w)). (18) 
Here AU (k) = f/(k, w) - £/ (k) . The first sum on the R.H.S. in Eq. (fill) can be estimated as 

i 2 k . 1 /n\V2 



1 /" d 2 k . 1 /n 

m / TTT n ° k + P - "sTJ ~ 

27lT 7|k|» Pj , l k l 7T 3/2 V ^ 



(19) 



(2tt) 7|k|> Pj 

Let us turn to the second term on the R.H.S. in Eq. (|18p . Since we are interested in H" 2 a {0,Pf) , let us substitute 
[G° CT (p + k, e + lj) + G arJ (p k, e — w)] « — u2+ ^2 ^ Pf) into this term. Then, one obtains the estimate for 
the second sum in Eq. ([TBI) 

- /*, / ^ <¥) !:HH - -a„v». ,20) 



(27r) J Jfc>p 



fc w 2 + (fc 2 /2) a 



This estimate results from the following hint. Let us introduce new variables uj/n 2 ^ 3 and k = k/n 1 / 3 . The integrand is 
maximal in the vicinity of the region uj/n 2 ^ 3 ~ k/n 1 / 3 ~ 1. Due to Eq. (|14p . one can put the lower limit equal to zero 
in integral (|20|) . Then the integral depends on the concentration as n 1 / 3 . The proportionality constant is calculated 
numerically, C 0.615. 



III. EQUATION OF STATE 

Taking into account condition ([T4"l) , one can neglect the contribution of expressions 1X0]) (jTTJ) , (fl~9| into the chemical 
potential as compared with ((20]). Thus, expressions (fTTj) and (|20|) alone contribute to the chemical potential. Therefore, 



H w 27Wn- Cn 1/3 . 



(21) 
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The chemical potential /z, given by Eq. ([2"Tj) . is a result of the first iteration. Let us now substitute expression (f2Tj) 
for the chemical potential into Eqs. © - (|5|) instead of /io and repeat all the calculations starting from Eq. (fTTj) . It 
is easy to find that, whether fj,Q or fi is taken as a starting point, the correlation contribution to the self-energy is of 
the same order of the magnitude given by (|20p . For this reason, the expression for the chemical potential (|2ip is the 
stable solution of the self-consistent system of equations, the constant C being of the order of unity. This chemical 
potential corresponds to the energy per unit volume 

e = nln 2 - -Cn 4 ' 3 . (22) 
4 

Then, at zero temperature the pressure is 

P = fin-e = nln 2 - j« 4/3 - (23) 
At small density n < (C/6irl) 3 ^ 2 , the compressibility dP/dV > 0, and the system is unstable. For higher density, 

3/2 3/2 

in the interval (C/6irl) < n < (C/AttI) , one has dP/dV < 0. The system is stable and has a negative pressure. 
The liquid state is reached at the upper limit of this interval 

n eq = (C/4irlf 2 . (24) 

where the pressure P = 0. 

Thus, for the electron- hole liquid in the equilibrium, the bound energy per one particle is 

C-3/2 1 

£ eq = -£ (n eq ) /n eq = (25) 

At the same time, for the isolated exciton the bound energy e ex is known to behave as s ex ~ —l/l. Since the excitons 
in the CQW experience repulsion,, the bound energy per particle is e ex > —l/l. Comparing this energy with e eq 
(see expression (|25p) one concludes that the electron- hole state possesses a smaller energy as compared to the exciton 
state for / 3> 1. In addition, note that the exciton in- plane radius R ex ~ l 3 ^ 4 and, thus, for the equilibrium carrier 
concentration obtained, one has n eq R 2 x ~ 1. Under such conditions, the exciton concept seems us to have no sense. 



IV. DISCUSSION AND CONCLUSION 



The results above are obtained for the case i > 1, i.e., I 3> ob- In this case, the ground-state of the system is the 
electron- hole liquid which possesses the equilibrium density ((24)) . This is a stable state which has the energy smaller 
than that of the exciton state. 

Consider, however, the case I <C 1. It follows from Eqs. (fTU|) . (fTTj) . ([2"U| that the kinetic energy Skin, the direct 
Coulomb interaction e e h and the intra-layer correlation energy e corr can be estimated as 

n 1/3 

V 

respectively. The equilibrium is reached at the concentration 



l + vl 



3/2 

(26) 



for which the total energy e eq is minimal and 



1/2 -j,if^/>l (97 , 

l + vl) ~\ -^ 1 /2 jif ^ <<1 • 

For the case I <C 1 considered in this place, the exciton bound energy e ex ~ — 1 in the units taken in the paper. 
(In fact, this energy is even larger due to repulsion between the excitons). Throughout the paper it is assumed that 
v 3> l.Thus, e eq e ex , as well as above in the case I 3> L Also, for the concentration (|26l) . the average distance 

between the carriers is of the order of n eq ^ 2 ~ {^^f^ ^ 1- At the same time, for I 1 the exciton in-plane 
radius R ex ~ Ob = !• Thus, the exciton in-plane radius i? e2; exceeds the average distance between the carriers, what 
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has no sense. Thus, both in the case I 1 and in the case I <C 1 the ground state of the system of spatially separated 
electrons and holes in the CQW is the electron-hole liquid. 

— 1/2 

As the intermediate case I ~ 1 is concerned, both e eq ~ — 1 and e ex ~ — 1. However, in this case n e q ~ Rex ~ 
as ~ 1 and, thus, the exciton concept seems has no sense as well. 

Let us return to the case I ^> l.An interesting feature of the results obtained is that both the bound energy ([24]) 
and the equilibrium concentration (|25[) of the electron-hole liquid in the CQW do not depend on the parameter v. 
This is in spite of the fact that formally this parameter was supposed to be large. This circumstance allows us to 
expect that the liquid equilibrium electron- hole state in the CQW may arise even for v ~ 1. However, for v ^> 1, 
the kinetic energy, which is of the order of ep ~ i, plays an insignificant role in the formation of the equilibrium 
state. The equilibrium within the electron-hole system is reached due to balance between the positive direct Coulomb 
interaction (1111) and the negative in-layer correlation energy (|20p . 

Let us pay attention to the interesting feature of the origin for the correlation energy associated with the RPA 
diagram. The correlation energy contribution results from the momentum transfer p ~ n 1 / 3 3> pf- This conclusion 
strongly differs from that of Cell-Mann and Brueckner who considered the many-body correlations for a 3D electron 
gas. Within their approach, the correlation energy originates from the momentum transfer p <C pp. Let us also point 
out additional arguments which approve the RPA approximation. While the Cell-Mann and Brueckner approach is 
justified for h>1, our approach holds for an arbitrary concentration n provided that v is large. 

Let us investigate how the results obtained can be associated with the experiments. Below, we follow the scenario 
proposed by Keldysh for the electron-hole liquid formation in semiconductors^. Consider a non-equilibrium inhomo- 
geneous electron-hole system in the CQW, composed of two coexisting phase: a rare exciton gas phase and a dense 
electron-hole liquid drops. Let R be the characteristic radius of the drop and r be the carrier (electron or hole) 
lifetime within this drop. Then, the loss of the carriers per unit time for the drop is — 7ri? 2 n eq /r. On the other hand, 
there exists an exciton flow j in the gas phase, induced by an external source. This flow gains the number of carriers 
in the drop per unit time +2TtRj. In the dynamic equilibrium, the loss and the gain balance each other, resulting in 
the dependence 

The concept of the electron- hole drop in the CQW fails if R < as = 1. Therefore, there exists a threshold value for 
the flow 

j>j c = ^, (28) 

which provides the formation of the drop. Thus, if the flow j is weak, the drop has no time to form. 

Such a scenario may account for the appearance of luminescent rings in the CQW in the series of experiments^—. 
These rings, sometimes, are associated with the exciton Bose condensate. However, in our opinion, they should be 
connected with the electron-hole liquid drop described above. First of all for these experiments, the criterion of the 
exciton existence R\ x n <C 1 seems to be violated. However, suppose that the core of the ring under the experiments 
is illuminated by the external optical pump resulting in the creation of the exciton gas flow j within the core. Then, 
the ring corresponds to the electron-hole liquid phase described above. It follows from Eq. (|28l) that the electron-hole 
liquid ring can dwell in the stationary state if the flow j exceeds a certain critical value. In the experiments, the 
luminescent ring appears if the optical pumping exceeds a certain power. Thus, the results obtained may be the basis 
for the explanation of these experiments. 

Let us mention that the correlation energy calculated above does not involve the effect of the exciton-like electron- 
hole correlation associated with the ladder diagrams in the electron-hole scattering channel. The involvement of this 
effect would result in the appearance of the nonvanishing anomalous exciton-like averages responsible for the creation 
of the dielectric gap on the Fermi surface^. However, the contribution of these diagrams into correlation energy is 
small in the parameters l/i/.This makes insignificant their influence on the thermodynamics of the system, the main 
subject of the paper. 

We are grateful to Yuri Lozovik for valuable discussions. 
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